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Abstract 

We extend the adiabatic bond charge model, originally developed for group 

IV semiconductors and III-V compounds, to study phonons in more ionic II- 

VI compounds with a zincblende structure. Phonon spectra, density of states 

and specific heats are calculated for six II- VI compounds and compared with 

both experimental data and the results of other models. We show that the 6- 

parameter bond charge model gives a good description of the lattice dynamics 

of these materials. We also discuss trends in the parameters with respect to 

the ionicity and metallicity of these compounds. 
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I. INTRODUCTION 



The adiabatic bond charge model (BCM) has been quite successful in explaining the 
phonon dispersion curves of group IV elemental semiconductorsEI and partially ionic III-V 
semiconducting materials! with a zincblende structure. In recent years it has been suc- 
cessfully applied to study phonons in semiconducting superlatticeJli, optical properties of 
open semiconductor surfaced^ and even sp^-bonded materials like graphite^ 
and fullerenes.i Quite recently, a modified version of the bond charge model has been appliedl^ 
to study the second order Raman spectra of AlAs and AlSb. 

In view of the similar dispersion curves of tetrahedrally connected III-V and II- VI ma- 
terials, it is surprising that no attempt to extend BCM to the latter has appeared in the 
literature. This may be partly due to the comments^ indicating that early attempts in 
this direction were not successful because it was foundS that in the case of II- VI materials 
the asymmetry of the bond charge position became too large to find a stable equilibrium 
position for the bond charges. These conclusions were based on the studies of the valence 
electron charge density using local psendopotentialil, which suggested a nearly complete 
charge transfer from the cation to the anion. These calculations actually overestimated 
the ionicity of these compounds and produced valence band spectra in strong disagreement 
with experimental photoemission results.^ Later, more accurate calculations using nonlocal 
pseudopotentialsQ showed better agreement with the experiments and yielded charge den- 
sities indicating a strong shift of the bond maximum rather than complete charge transfer. 
In fact, the charge density plots for III-V and II- VI compoundslll are nearly identical except 
that the charge maxima in the latter appear to be slightly shifted toward the anion, indicat- 
ing that, despite their greater ionicity, II- VI compounds are dominantly covalent in nature. 
In view of these results, it is expected that the BCM should give a good account of the 
phonons in II- VI materials provided the bond charges are placed at physically reasonable 
places suggested by the pseudopotential calculations. 

Traditionally, the lattice dynamics of these materials has been done using rigid ion or 
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shell models.cJllj These models give good fits to the observed phonon dispersion curves at 
the cost of a large number of adjustable parameters (10 or more), some of which have no 
physical interpretation. Recently, ab initio calculations of phonon spectra have appeared,^ 
but they are are not feasible for studying large systems such as alloys like Cd^^Hgi-ajTe or 
thick superlattices. Therefore it is desirable to have a realistic model with fewer, physically 
meaningful, parameters that is easy to extend to more complex systems. In this paper 
we show that the 6-parameter BCM provides a good description of the phonons and other 
lattice dynamical quantities such as elastic constants and specific heat in II- VI materials. 

The rest of the paper is organized as follows. In Sect. II we provide a brief overview of 
the bond charge model.0^ In Sect. Ill we discuss the results for six II- VI compounds and 
in Sect. IV trends in the parameters are discussed. 



II. ADIABATIC BOND CHARGE MODEL 

The adiabatic bond charge model (BCM) for homopolar semiconductors^ and partially 
ionic III-V compoundJi is the simplest empirical lattice dynamical model that correctly 
describes the phonon dispersion curves of covalent crystals. In the BCM the valence electron 
charge density is represented by massless point particles, the bond charges (BC's), that follow 
the ionic motion adiabatically. The BCM unit cell consists of two ions and four bond charges 
that are placed along the bonds between the ions. In homopolar covalent crystals the bond 
charges are placed midway between the neighboring atoms while in III-V compounds the BC 
divides the bond length in the ratio of 5:3. This is consistent with nonlocal pseudopotential 
calculations for the valence electron charge density0 that indicate that the charge density 
maximum in III-V compounds shifts toward the group V element. This shift is even stronger 
in the case of II- VI compounds, reflecting their more ionic character. The BCM parameter 
p which measures the polarity of the bond is defined in terms of the ratio in which the 
BC position divides the bond length. If t is the bond length and ri = (1 + p)t/2 and 
r2 = (1 —p)t/2 are the two ion-BC distances then p = (ri/r2 = 1) for homopolar materials 
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and p = 0.25 (ri/r2 = 5/3) for III-V compounds. In our extension of the BCM to II- VI 
materials we have chosen to use p = 1/3 corresponding to the ratio ri/r2 = 2 which is based 
on the results^ of microscopic calculations. This choice will be discussed in more detail in 

Sect. |rv[ 

The cation and the anion interact with one another and with the bond charges via central 
potentials (puit), 0i(ri), and 02(^2), respectively. The bond charges centered on a common 
ion interact via a three body Keating potential,^ V^f, = B„(X^-XJ + a^Y /Sa^, where is 
the distance vector between ion cr (a = 1, 2) and BC i, B^r is the force constant and is the 
equilibrium value of ■ XJ\. The bond charges centered on a particular ion also interact 
directly with one another through a central potential, where r^^^ is the distance 

between the bond charges centered on the cation {a = 1) or the anion (cr = 2). Finally, the 
ions and the EC's interact via the Coulomb interaction characterized by a single parameter 
Z'^/e where —Ze is the charge of a BC, and e is the dielectric constant. Each of the ions is 
presumed to have a charge +2Ze so that the net charge in the unit cell is zero. 

To reduce the number of parameters it is assumed^ that ip[ = 'ip'2 = 0, ip" = —ijj!^ = 
{B2 — Bi)/8 and (1 + p)(f)'i + (1 — p)02 — 0- If ^se these constraints on the total lattice 
energy per unit cell 

$ = A[Ut) + 01 (n) + Mr2)] - + 6[n6 + + Mrl^) + V^2(rJ?)] (1) 

along with the equilibrium conditions, d^/dt = and d^/dp = 0, we findS 

J,' ^ c 

0n = — ^ 

01 ^ r, daM l-p Z^e'^ 
Vi dp 1 + p e 

r2 dp 1 — p e 

The conditions for stable equilibrium, d^^/dt^ > and d^^/dp^ > 0, further yield 

0'.' , ,2.0'/ B2, , ,2.02 128 ^2 

4f + (1 + + f ) + (1 - + f ) - 5^a.„- > (3) 

and 



The Madelung constant of the model is defined by writing the Coulomb energy per unit 
cell as — (y.Mi^Ze)^ / €t. For p — 1/3 the values of Q-mi dcum/dp and d'^o^M/dp'^ are found 
numericallyii to be 5.0598, 4.0539 and 17.46, respectively. In Eq. (^) and Eq. (^) the force 
constants are in units of e^/va, where Va is the unit cell volume. 

With 'ip'i, 4^2^ i^i and ip!^ given as above, the six free parameters of the model 

are (p"^, 0'/, Bi, B2 and Z'^/e, which we adjust to fit the neutron scattering data and 
the measured elastic constants. The phonon eigenfrequencies and eigenvectors are found by 
diagonalizing the dynamical matrixS constructed from the BCM equations of motion 

Mc.;^^ = [R + 4^^C«]u + [T - 2^^C^]v 

= [T+ - 2^Cj]u + [S + ^C,]v (5) 

Here M is the mass matrix for the ions and u and v are the vectors formed by the 
displacements of the ions and the EC's, respectively. The matrices R, T, and S are the 
dynamical matrices for the short range ion-ion, ion-BC and BC-BC interactions and C^, 
and Cg are the corresponding Coulomb matrices which are evaluated by Ewald's method.!! 
Explicit forms of R, T, and S can be found in the appendices of Refs. ffl and H. 



III. RESULTS AND DISCUSSION 

Figure 1 shows the dispersion curves for six II- VI materials along with the existing 
neutron scattering data. Figure 2 shows the corresponding densities of states (DOS). We 
present in Table I the BCM parameters and in Table II the calculated and measured elastic 
constants. In all cases the overall agreement with the experimental data is fairly good and 
is of the same quality as for BCM fits for III-V compounds.il For comparison we have also 
included the dispersion curves for GaAs and InSb calculated using the parameters from 
Ref. 1^. For CdTe our six parameter fit is as good as the 14 parameter shell model fit of 
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Rowe et al^ and the 11 parameter rigid ion model fit of Talwar et al^ All three models 
show a slight upward bend in the TO branch in the (100) direction, which is in contrast 
with the results of a recent ab initio calculation.El 

The three models give very similar predictions for the phonon density of states. The 
principal difference is that the shell model does not predict a gap in the DOS between the 
acoustic and optical contributions, while the BCM has a smaller gap than that seen in the 
rigid ion model. This underscores the fact that the BCM is, in many ways, an intermediate 
model between the shell model and the rigid ion model. The shell model takes care of the 
electronic polarizability explicitly by attaching deformable shells to the ions. The BCM 
partially accounts for the electronic polarizability through the adiabatic motion of the bond 
charges, while the rigid ion model ignores it completely. 

For ZnS and ZnTe our fits are comparable with the 10-parameter Valence Shell Model 
results of Vagelatos et al^ For ZnTe the BCM predicts a large dispersion in the LO branch 
near the zone edge. However, the maximum deviation from the experimentally measured 
frequency at the X point is only about 7 percent. For both ZnS and ZnTe the shape of the 
optical branches in the (110) direction is different from the results of Ref. ^ but is similar 
to that predicted^! by ab initio calculations. For ZnSe also the agreement with the neutron 
data@ and the measured elastic constants^ is fairly good. 

Mercury compounds, because of their semi-metallic nature, deserve a separate discus- 
sion. Because of the zero band-gap, the energy for electronic transitions from the valence 
band to the conduction band is comparable to the optical phonon energy. For HgTe Raman 
measurements@ at 90 K and infrared reflectivity measurementsil at 77 K yield cj^q ~ 138 
cm~^ whereas infrared spectra at 8 K gaveil cj^o ^ 132 cm~^ . This difference was at- 
tributed to the large number of carriers at higher temperatures. The same reason is invoked 
to explain the degenerate values of u^q and Uj^q found in neutron scattering experiments.^ 
Because of this controversy we did not use optical phonon frequencies near the zone center 
in our flt for HgTe. It is seen that the agreement with the acoustic and transverse optical 
phonons is good. However, the flt for the LO branch is not of the same quality, although 



the deviation from the experimental points is only a few percent. In the (111) direction 
the BCM predicts that the LO branch will dip downward instead of the upward trend ob- 
served experimentally.!^ For HgSe the agreement with the available neutron data on acoustic 
phononJ^ and optical measurement^ and measured elastic constants^ is very good. Be- 
cause of the lack of neutron data on optical phonons we cannot comment on the accuracy of 
the optical branches. However, it should be mentioned that the BCM and the 11 parameter 
rigid ion modefl give similar behavior for the optical branches. 

To further check the parameters we have calculated the specific heats@ for all the six 
materials using the parameters given in Table I. The results are shown in Fig. 3 as plots of 
log C vs. log T along with the experimental data. In every case good agreement is obtained 
with the experiments, giving further support for the parameters used and the calculated 
density of states. 

We have thus demonstrated that the 6-parameter BCM provides a good description of 
the phonons and other lattice dynamical quantities such as elastic constants and specific 
heat for II- VI compounds with zincblende coordination. The overall agreement with the 
neutron data is very good. However, some discrepancies remain, particularly in the LO 
branch near the zone edge, where the BCM predicts a large dispersion in almost every 
material including the III-V compounds. This is most obvious in HgTe and is indicative 
of the failure of the BCM to account for the polarizability of the ion core. Because of the 
associated macroscopic field, the LO phonons are more affected than the other branches. The 
calculatedill static dielectric function e{q) for III-V and II- VI compounds is known to have 
considerable structure, so including a charge form factor should remedy this discrepancy.0 



IV. TRENDS IN PARAMETERS 

Some trends in the parameters presented in Table I are immediately obvious. One 
notices that as one goes from group IV elements to III-V compounds to II- VI compounds 
the parameters involving BC's change uniformly. For group IV elements the bond charge is 
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situated midway along the bond and the ion-BC and BC-ion-BC force constants are equal 
for the two ions. However, in III-V compounds the BC shifts toward the anion which results 
in higher values for 02/3 and B2 than 0'//3 and Bi respectively. This trend continues as 
we move to II- VI compounds in which the BC is even closer to the anion. This pattern in 
the values of these parameters can be traced to the equations linking (f)'^ and 02 to and 
daM/dp. It should be noted that the values for and dau/dp for II- VI compounds are 
higher than those for III-V compounds. Apart from these obvious features, there are no 
other discernible trends in the parameters with respect to the ionicity or the bond length. 
However, it is seen that the ion-bond parameters 0i/3 and 02/3 are considerably lower 
for mercury compounds than for other materials. This is reasonable in view of the semi- 
metallic nature of these materials and the fact that these parameters represent off-diagonal 
contributions to the dielectric function. 

The effect of ionicity on the phonon dispersion curves can be investigated by studying 
the isoelectronic sequence of materials in which the bond lengths and the average mass in 
the unit cell are almost same. Two such sequences are Ge-GaAs-ZnSe and aSn-InSb-CdTe. 
Increased ionicity results in a general lowering of all frequencies and elastic constants and a 
hfting of the LO-TO degeneracy at the V point and the LO-LA degeneracy at the X point. A 
glance at the BCM parameters for these materials shows that the only pattern is a decrease 
in the magnitude of the bond charge Z and, in the case of Ge-GaAs-ZnSe, a decrease in 
0-'j/3 with increased ionicity. For the other sequence, 0i'j/3 is almost same for a-Sn and 
InSb but it is smaller for CdTe. A cross comparison of the corresponding materials in the 
two sequences shows that moving down the periodic table, with its concomitant increase in 
metallicity, yields an increase in the ion-ion interaction 0fj/3 while the parameters involving 
the bond charges decrease. 

We should also comment on the choice of the equilibrium positions of the bond charges. 
In principle, p should be treated as the seventh adjustable parameter of the model. However, 
we decided to use the physically reasonable value of 1/3 for p. This corresponds to dividing 
the bond length in a ratio 2:1 and is consistent with the pseudo-potential calculations of 
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the valence electron charge density.EJ However, we were also able to find values for the six 
parameters which still satisfied the stability conditions (0) and (^) and which gave satisfac- 
tory fits for p as high as 0.55. The parameters that varied most with p were the ion-BC 
force constants; increasing p led to a larger 02/3 and a smaller 0'//3. The remaining pa- 
rameters also changed, but by much less. In every case, the acoustic phonon curves showed 
very good agreement with the neutron scattering data. However, the agreement with the 
optical phonons slightly worsened as p increased. These results highlight the arbitrariness 
involved in defining the equilibrium position for the bond charges. Our choice of p = 1/3, 
which coincides roughly with the position predicted by the pseudopotential calculations^, 
still gave the best overall agreement with the experimental results. 



V. SUMMARY AND CONCLUSION 

We have applied the adiabatic bond charge model to study phonons in six II- VI com- 
pounds with a zincblende structure. The theoretical predictions of the 6-parameter BCM 
are in good agreement with the available neutron data and the experimentally measured 
elastic constants and specific heats. Some minor discrepancies in the LO branch near the 
zone edge are believed to be due to the incomplete description of the electronic polarizability 
of the ions. These deviations are larger, though still only a few percent, in the case of HgTe 
as expected from its semi-metallic nature and consequently stronger screening effects. In 
conclusion, we have found that the 6-parameter adiabatic bond charge model provides a 
satisfactory description of the lattice dynamics of tetrahedrally connected II- VI compounds. 
The agreement with the experimental data is of the same quality as for III-V compounds. 

We have also discussed some broad trends seen in the parameters. We find that the 
parameters for the bond charge-cation short range interactions decrease with a corresponding 
increase in the parameters for the bond charge-anion interactions as one goes from group 
IV elements to III-V to II- VI compounds. 

We thank T. Golding for useful conversations. This work was supported by the National 
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FIGURES 

FIG. 1. Calculated phonon dispersion curves for CdTe, ZnS, ZnTe, ZnSe, HgTe, HgSe, GaAs, 
and InSb. The BCM Parameters for GaAs and InSb were taken from Ref. |3|. Empty circles indicate 



neutron scattering data taken from Refs. (CdTe), (ZnS and ZnTe), (ZnSe), 28 (HgSe and 



HgTe), M (GaAs), and 37 (InSb). 



FIG. 2. Phonon density of states for II- VI compounds calculated using the root sampling 
method. The fine structure on the curves is an artifact of the numerical method. 

FIG. 3. Log C vs log T plots for the calculated and measured specific heats of several II- VI 
materials. The experimental data is taken from Refs. p^41. 
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TABLES 

TABLE L BCM parameters for group IV elements, IILV and ILVI compounds. Force constants 
are in units of e^/fa, where Va is the unit-cell volume. 







0'iV3 


<A^73 


Bi 


B2 








6.21 


6.47 


6.47 


8.60 


8.60 


0.1800 


1.47 




6.61 


5.71 


5.71 


8.40 


8.40 


0.1620 


1.61 


a-Sn*^ 


7.43 


5.59 


5.59 


7.80 


7.800 


0.163 


1.98 


AlAs^ 


5.80 


2.27 


15.48 


5.79 


8.54 


0.1800 


1.21 




6 04 


2.4 


1 7 91 


5 20 


1 


2030 


1 36 


GaAs^ 


6.16 


2.36 


16.05 


5.36 


8.24 


0.1870 


1.43 


GaSb^ 


6.77 


2.37 


13.10 


6.28 


7.08 


0.1600 


1.52 


InP^ 


7.16 


2.95 


21.62 


3.43 


8.37 


0.2490 


1.55 


InAs*" 


7.31 


2.64 


17.86 


3.99 


7.30 


0.2100 


1.60 


InSb^ 


7.47 


2.33 


14.09 


4.56 


6.24 


0.1720 


1.64 


ZnS 


5.74 


0.79 


29.90 


0.83 


15.40 


0.2130 


1.05 


ZnSe 


5.01 


1.19 


22.82 


1.21 


15.65 


0.1790 


1.03 


ZnTe 


5.51 


1.06 


22.93 


1.07 


17.00 


0.1800 


1.05 


CdTe 


6.85 


0.77 


23.34 


0.39 


15.44 


0.1830 


1.15 


HgSe 


5.32 


0.15 


14.01 


0.35 


17.50 


0.1095 


0.91 


HgTe 


6.46 


0.081 


13.46 


1.08 


15.60 


0.1062 


1.03 



^Parameters from Ref. |^. 
^Parameters from Ref. |^ 

■^eoo for HgSe from Ref. |32l for HgTe from Ref. 27 and for the rest of the materials from Ref. 
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TABLE II. Theoretical and measured values (in parantheses) for the elastic constants in units 
of lO"*^^ dynes/cm^. 



Cll 



Cl2 



C44 



ZnS^ 
ZnSe^ 
ZnTe^ 
CdTe'^ 
HgSe= 
HgTe'^ 



10.907 (10.46) 
8.996 (8.59) 
7.138 (7.13) 
5.675 (5.35) 
6.218 (6.22) 
5.631 (5.631) 



'^Measured values from Ref. |3^. 
'^Measured values from Ref. B4. 



6.498 (6.53) 
5.064 (5.06) 
4.233 (4.07) 
4.073 (3.68) 
4.647 (4.64) 
3.785 (3.66) 



4.678 (4.61) 
4.056 (4.06) 

3.122 (3.12) 
2.047 (1.994) 
2.262 (2.27) 

2.123 (2.123) 



"Measured values from Ref. 



"^Measured values from Ref. 35 
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